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Boost modes for a massive fermion field 
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We have shown that Wightman function of a free quantum field generates any complete set of 
solutions of relativistic wave equations. Using this approach we have constructed the complete set 
of solutions to 2d Dirac equation consisting of eigenfunctions of the generator of Lorentz rotations 
(boost operator) . It is shown that at the surface of the light cone the boost modes for a fermion 
field contain 5- function of a complex argument. Due to the presence of such singularity exclusion 
even of a single mode with an arbitrary value of the boost quantum number makes the set of boost 
modes incomplete. 
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I. INTRODUCTION 
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> 

O, . 

^^ , It is well known that quantization procedure in a quantum field theory implies expansion of the field operator in 

terms of a complete set of modes which are solutions of the corresponding (Dirac or Klein-Fock-Gordon) classical field 
'• ' equation. Therefore finding exact solutions for these equations is one of the central points of the quantum field theory 
in the presence of a classical background. 

One of the most powerful instruments for solving partial differential equations is based on symmetry properties of 
the physical system described by the equation. According to the Noether's theorem [1] any differentiable symmetry 
of the action of a physical system has a corresponding conservation law. One can always choose a set of variables in 
^H such a way that a change in one of them corresponds to the symmetry transformation. Since the generator of the 
~^ ' symmetry transformation commutes with the differential operator of the equation, the variable corresponding to it 
can be separated. If there are enough symmetries for the equation, such that their generators commute with each 
other, the problem of finding the solution can be reduced to solving, generally speaking, of a second order ordinary 
^N ' equation. The set of such solutions is labeled by the eigenvalues of the generators and is complete. 

All most important and widely used solutions for relativistic quantum equations in the presence of classical ex- 
ternal fields were obtained using this method. We mean the Coulomb field, see, e.g., [3, the constant magnetic [3| 
and electric [3,0 fields, the field of a plane electromagnetic wave [g|, and some other fields of more sophisticated 
f^ configurations [J, 0] . In all these cases the solutions were eigenfunctions of linear combinations of generators of time 

1^ * _ and spatial translations or spatial rotations, and were labeled by values of components of either 4-momentum, or 

f-^ ' angular momentum respectively. However, the Poincare group, the group of isometries of Minkowski spacetime (MS), 
includes also Lorentz rotations (boosts). The boost symmetry was almost never used for field quantization. Certainly, 
this is because the boost generator does not commute with the Hamiltonian, and thus the boost quantum number 
and energy cannot be parts of any complete set of observables simultaneously. Nevertheless, the boost modes can be 
very useful, especially for the quantum field theory in a curved space, or if the symmetry of MS with respect to time 
and/or space translations is broken in the presence of a classical background. In these cases the boost symmetry can 
5^ ] appear to be the sole symmetry for the quantum field. Dilaton gravity in two dimensions [3] , as well as Schwarzschild 
geometry Q, are the examples. Therefore analysis of the properties of boost modes is important. 

For the first time the boost modes for a free massive scalar field were discussed in Ref. Q by W. Unruh, though 
their explicit form was not specified. Then some properties of scalar boost modes were studied in Refs. |10l - [l3 |. In 
particular, the remarkable properties of boost modes in 2d MS were ascertained in Ref. [l2|- It was shown that (i) 
the zero boost mode of a free massive scalar field coincides up to a trivial constant factor with the positive-frequency 
Whightman function, and (ii) the boost modes considered as functions of the boost quantum number possess the 
Dirac (5-function singularity at the surface of the light cone. The properties of the boost modes for the case of a 
free massless fermion field appeared to be even more interesting. It was shown in Ref. [13[ that boost modes of 
two-dimensional massless fermions on a light cone are expressed in terms of the delta function of a complex argument. 
Boost modes of a massive fermion field have been first considered in Refs. [ij, [la], see also Ref. 16]. However the 
results of Refs. |14|-[l6| contain some discrepancies which we discuss in Sec. HI of the present paper. 

In this paper we consider boost modes for a free massive fermion field. To construct them we start from the 
Wightman functions. The positive and negative frequency Wightman functions for free fields are explicitly determined 
(with accuracy to a constant factor) by Lorentz and translational invariance of the theory [l7[ • And vice versa, once 
the Wightman functions are known, one can use those symmetries to construct any complete set of positive and 
negative frequency solutions of Klein-Fock-Gordon (KFG) or Dirac equations. In fact, this is the direct consequence 



of the Wightman reconstruction theorem [T^ [l^ . 

Since the Lorentz rotation singles out a two-dimensional plane in MS, we will discuss the specific properties of boost 
modes by examples of free 2d KFG and Dirac equations. In the next section we will realize the above formulated 
approach by the example of a single-component massive neutral field. It is worth noting that, as opposed to the case 
of plane wave modes, analytic properties of boost modes change dramatically when one passes to multi-component 
fields which are considered in Sec. III. The discussion of the results and conclusions are given in Sec. IV. 

II. BOSONS 

In this section we will consider the case of a free neutral massive scalar field in two-dimensional MS. We will start 
from two-point Wightman function [l8| which for a free field theory coincides with the positive-frequency part of the 
commutator of two field operators (Pauli- Jordan function in 4d theory of a scalar field) . Positive- frequency Wightman 
function A''"''^(a;), x — (i, z), for a massive scalar field satisfies KFG equation^ 

/CA(+)(x)=0, IC = d^-dl+m^, (1) 

contains only positive frequencies and is invariant with respect to Lorentz rotations (boosts) 

6;c A^+' (x) - 0, 6yc = t{zdt + td,). (2) 

These conditions determine A(+)(a;) accurate within a constant factor. Indeed, any positive-frequency solution of 
KFG equation can be written as follows 

$(+)(x)= f d^p(j>ip)S{p^ ^m^)9{p'^)e-'P'=, (3) 



where <j>{p) is a certain function of 2-vector p = (p°,p^), 6'(p°) is the Heaviside step function. After the change of 
variables 

p = ^coshq, p = fisinhq, 

where q is rapidity, and integration over ^ we obtain the equation 

CO 

$(+)(a;) = i /" dg0(m cosh q,msinhq)e~™(*'=°''''«-^ '''''''«). (4) 

— oo 

It can be easily seen that condition ([2]) is satisfied only if <j){q) — const. Choosing (p{q) ~ i/27r we arrive to the 
standard representation for A(+'(a;), compare, e.g., [H, [12 



I 



A(+)(a;)^_ ^q^-im{t cosh q^zsinhq) ,^^ 

It is assumed in ([S|) that an infinitely small negative imaginary part is added to t [12]. Another representation for 
A(+)(x) reads, see, e.g., [13], 

OO 

* f dp __, 



A(+)(a;) = — / JLe-'e,t+^pz ^ ep = ^p^+m\ (6) 

47r 

from now on we shall omit the index of the spatial component of the 2-vector p. 

According to the Wightman reconstruction theorem [l8| the two-point Wightman function uniquely determines the 
quantum theory of a free field. Particularly, it allows to reconstruct any orthonormalized and complete set of solutions 



^ The natural units h = c = 1 are used throughout this paper. 



for the field equation. Indeed, it follows from the translational invariance of the theory that Wightman function with 
an arbitrarily shifted argument 

A'^+\x-u), m = {m°,u1} (7) 

satisfies KFG equation. Functions (O constitute an overcomplete set of modes since they are labeled by two indepen- 
dent parameters u'^,w^. Anyway, any complete set of positive-frequency solutions Fa{x) 



daFa{x')F:{x") = -iA^+\x' - x"), (8) 

orthonormalized by the condition 

{Fa,Fa-) = i I F:{x)dtFa.{x)dz^5{a-a') (9) 

can be represented as 

Faix) = / d''ufa{u)A''+\x - u) . (10) 



However, since the set of solutions (O is overcomplete, the coefficient functions fa{u) cannot be determined uniquely. 
To make the choice of fa (u) unique one should impose a restriction on u which can be chosen for reasons of symmetry. 
First, we will illustrate this procedure for the trivial case of plane waves Op{x) which are eigenfunctions of the 
generator of spatial translations 

-id,ep{x)^pQp{x), (11) 

In this case it is convenient to limit the family of A*^"*'-' {x — u) only by functions spatially shifted with respect to each 
other. This means that coefficient functions fp{u) have the form fp{u) = d{u'^)'dp{u^), so that 

oo 

ep{x)= [ du^dp{u^)A^+\t,z-u^). (12) 



Let us substitute this expansion into Eq. (TTlT) . Taking into account that Wightman function A'^+'(i, z) exponentially 
tends to zero when \z\ — > oo, we can integrate the left-hand side of the resulting equation by part and arrive to the 
following equation for -dpiu^) 

-id'p{u')=pi)p{u'), 

so that dpiu^) = const • e*''" . Substituting now dp{u^) of that form into ([T2|) and using ([6|) we finally obtain after 
integration over u^ 

ep(x) = ^^e-'^'''+'p^ (13) 

V47rep 

The normalization constant here was determined from the condition ([9|). 
Functions Op{x) satisfy the relation 

oo 

/ dpep{x')Q;{x") = -iA^+\x' - x"), (14) 

— oo 

and thus the set P^ is complete. 

We are interested in the set of boost modes ^k ^(x) which are eigenfunctions of boost operator B/c 

B,c^i+Hx) ^ ^¥+Hx) . (15) 



It is worth noting from the very beginning that Wightman function A(+)(x) is a zero mode ^q '{x) of this set as it 
follows from Eqs. ^^ and ©. 

To obtain the boost modes we will confine u in ([TU|) to one of the orbits of the restricted Lorentz group 

u^ — vP — u^ = ±w^ , V — const , (16) 

so that f{u) ^ 6{u'^ =F u^). Choosing the upper sign in (fT6|) for definiteness, we put 

M = liq = (w coshq, usinhq), (17) 

and rewrite Eq. (jlOp in the form 



^^(a:) = j dgCK('7)A(+)(x - u,). (18) 

— oo 

Substituting (TTSl) into Eq. (|T5|) and using the relation 



_9_ 



B^A<^+\x - u,) ^ ~^i—A^+\x - ug) , (19) 



which can be straightforwardly obtained with the help of representation ([SJ, we get the following equation for Ck('?) 

Hence, Ck('z) = const e~""', and we finally get from p^ with due regard for condition © the following representation 
for the positive frequency boost modes 

oo 
^(^+)(a;) = _1_ f ^q^-'^m{t cosh q-zsinhq)-ZKq^ ^20) 

— oo 

which was earlier obtained in Refs. |10l - [l3 | by another method. The negative frequency boost modes ^k (x) are 
defined [l^] as 

M'M(a;) = vl/L+i*(2;). (21) 

Boost modes (|20l) . ((2T|) are distributions with respect to both x = {t, z) and spectral parameter k. They are defined 
on the class of smooth functions of rapid enough descent. 

The modes ([201) 7 (EH) constitute a complete set since they satisfy the condition 

00 

/ dK^i^\x')-^'i^^*{x") = TiA(±)(a;' - x"), (22) 

—00 

and hence can be used as a basis for quantization of a neutral scalar field, see Ref. [12]. Since the sign of a particle 
energy is Lorentz invariant, the vacuum states in the boost and plane-wave quantization schemes are identical. Hence, 
these two quantization schemes are unitary equivalent, see Ref. |12| . 

The most remarkable property of boost modes is their behavior at the light cone. It is easily seen from ([20)1 . (PT|) 

that at the vertex of the light cone ^L (x) possess the Dirac (5-function singularity 

*W(0) = i=<5(^). (23) 

It was shown in Ref. [lj| that ^^ (x) possess a (5-function singularity at the lines x± = t ± z = as well. Therefore 
contribution of the single spectral point k = to physical quantities can be finite. We will illustrate this by the 
example of Wightman function. 



Taking into account the property of translational invariance we can rewrite Eq. (j22p in the form 

oo oo 

A(+)(x' - x") = i f dK.^';+\x')^i+^* {x") = i f dK.¥+Hx' - x")¥+^*{0) . (24) 

— oo — oo 

Using now Ea.p3| we obtain 

A(+) {x' -^x") = ^ *^+) {x' - x") , (25) 

i.e. the integral over k in (I24p is determined by the point k = entirely^. 

This resuh means that the point k — cannot be deleted from the spectrum, or in other words the integral over k 
in (P^ . (|24l) cannot be changed by its principal value 



dn... T^P.v. / dK... . (26) 

— oo 

Thereby, the family of boost modes does not constitute a complete set in MS after excluding the zero mode. However, 
the equality in Eq. (|^ could be restored if we cut the light cone out of MS, compare [l2|- This is because all the 
points where boost modes possess S{k) singularity are located just at the light cone. 
Furthermore, the Wightman function Eq. (|24p cannot be represented in the form 



A(+)(x' - x") = ^ / dK{^i+\x')^i+^*{x") + *L+i(a;')*L+i(x")*} • (27) 



If it were so, we could rewrite (j24p . due to translational invariance of the Wightman function and the boost mode 
property (|23p . as follows 



A(+)(a;' - x") ^A= dn{^'^+\x' - x") + *L+i(2;' - x")}5{k) . (28) 



But Eq. ([28|) is evidently meaningless. Indeed, the distribution 5{k) is defined on the class of functions continuous at 
the interval including the point k = 0. Therefore Eq. (P5| should be understood as 

oo 

A^+\x' - x") = A= I dK{-^'^+\x' - x") + *L+J(.t' ~ x")}e{K)5{K) , (29) 

— oo 

where 9{k) is the Heaviside step function. However, the product of two distributions 6{k)6{k) is not defined. 

The authors of Ref. [21j do not agree with this statement. Discussing their Eq. (2.129), which in our notations 
coincides with ([27]), they admit that this expression is undefined if either of the two points x' , x" is located on the 
light cone. They think that smearing of the distribution A'^^ [x! — x") with compactly supported functions f{x') and 
g{x") improve the situation and the smeared "Wightman function" (2.131) is well defined on the whole MS. Our 
analysis shows however that this is not correct. Indeed, if we use Eq. ([?7|) in the form (^5]) . we see that this expression 
is undefined for arbitrary x', x" . Moreover, smearing cannot improve the situation since it does not influence the 
(5- function in ([25]) . The correct expression for the smeared Wightman function can be easily obtained from Eq. (|24p . 
It reads 

oo 

A^+\f,g) = j d^x'd^x"r{x')g{x")^^+\x'~x")^t j dK^i+\f,g)^^+^\Q) = ^M>^+\f,g), (30) 



We have mentioned already that Eq. I|25| l directly follows from Eqs. {J} and ||15|I accurate within a constant factor. Then Eq. (|23| l could 
be derived from Eqs. II25I I. II24I I even not using the representation I I20I I. 



^, = 




a) 



b) c) 

FIG. 1: Classical trajectories of massive boost particles. 



*l+n/,.9) - / d^x'd^x"rix')gix")^(+Hx' - x") , 



(31) 



where 



is the smeared boost mode. We see that again only the single spectral point k = contributes to the integral (l30l) for 
the smeared Wightman function. 

Note, that we could try to use expression (|27|) directly, not applying the property of translational invariance to it. 
In that case we should attach exact mathematical meaning to (j27p first. Namely, we should write it in the form 



E-i-O 



A^+\x',x") = ilim / dn{^^+\x')^^+^*{x") + *L+i(a;')*Lt,^*(x")} • 



(32) 



But it can be easily seen that function A(+'(x',x") in (|32|) is not translationally invariant and thus has nothing to do 
with Wightman function for a free field in MS. Indeed, if it were translationally invariant, it could be written down 
as 



^^+\x',x") = ^ lim f dK{¥+Hx' - x") + *L+i(x' - x")}d{K) , 

V2 e->0 J 



(33) 



and would be equal to zero identically. 

Singular behavior of the boost mode at the point k = can be understood in terms of classical trajectories of free 
particles with a given value of the boost parameter k 



K = ZSp — tp, 



(34) 



where p = mzj \J\ — z^ is the momentum, Sp = m/\/\ — i^ - the energy of a particle, z = dz/dt. Let us rewrite 
Eq. ([311) in the the form 

(35) 



K I 

z = zt^ VI - z^- 

m 

After differentiating both sides of Eq. (|35l) with time one can easily find its regular 

z — vt -\ vl — z)-^, V — const , |w|<l, 

m 

and singular 



(36) 



z = sgn(K)Wt2 + 



(37) 



solutions. 

When K 7^ the singular solution is represented by two branches of hyperbola ([57]). The branch with k > is 
located in the right wedge of MS, while the branch with k < in the left one, see Figs. la,b. Free boost particles 
with the given value of k propagate along the regular world lines which are different straight lines tangent to the 
corresponding branches of hyperbola ((37)) . It is worth noting that classical particles with positive values of the boost 
parameter k cannot penetrate to the left wedge, as well as particles with k < cannot penetrate to the right one. 

The world lines of boost particles with k = are given by 

z = vt. (38) 

All of them cross the vertex of the light cone /iq and any world line crossing /ig belongs to the family of trajectories 
with K = , see Fig. Ic. Hence, elimination of the point k = from the spectrum is equivalent to prohibition for 
particles to cross the point /lo, or to pricking it out of MS. Besides, deleting the point k = from the spectrum, we 
delete a bunch of infinite number of trajectories, thus loosing a substantial part of degrees of freedom. This means 
that the point k = is of nonzero measure at /iq and explains the (5-function singularity of the boost mode there. 

The points lying at the light cone surface need a separate consideration. The lines z — zLt cannot be trajectories of 
massive particles. But the branches of hyperbola (1571) degenerate into lines constituting the surface of the light cone 
in our 2d problem when k — > 0. As it can easily be seen from Figs. la,b, all regular trajectories except those which are 
tangent to the hyperbola branches at the point t — also tend to the lines z = zLt in the limit k -^- 0. Hence, these 
lines z = ±t are the lines of condensation for regular trajectories with k — > 0. Possibly this explains the existence of 
(5-function singularity of boost modes at the light cone surface. 

It is worth noting that there exist another representation for the boost modes (|20p . (PT|). Keeping in mind that 
an infinitely small negative imaginary part —ia is added to t in (j20p . after the change of the variable of integration 
q = \nu we obtain for "^n [x) the following expression 

*«^^ (^) = W^ / rf"""^^"'e" " """ - (39) 



with /3 =^ — (a =b 2X+) and 7 = ^(cr ± ix_). Using now formula 3.471(9) of Ref. l22j we get 

-^i^Hx) = — -= — K,^{w±) , w± = m^e±*'^(x+T«cr)(a;_T«cr), (40) 

7rV2 \x+ T icrj 

where Kii^{w±) is the Macdonald function, and distributions (i^ =F ia)^ should be understood as 

(e T ^<J)^ = ^0(0 + e^"\'O^0{-O, (41) 

see Ref. [H]. Using Eq. (gl]) we obtain 

^.{x) = ei-x+)ei-x^)¥p{x) + eix+)ei-x^)¥,''Hx) + e{x+)e{x^)¥f^x) + ei-x+)eix^)¥^^Hx), (42) 

see [lll[i2|i where 



(43) 

are the expressions for the boost modes in the past (P), right (i?), future (i^) and left (L) wedges of MS respectively. 
Here H- ' (w) are Hankel functions. 

Putting K = in Eq. (|43p and taking into account Eq. (|25|) one can easily reproduce the Wightman result [20| for 
A'+^(x). Applying the same procedure to Eq. (PII|) we get 

A^^Hx) ^ ±^Ko{w±) . (44) 

This is a new compact representation for the Wightman function in 2d scalar theory. 



8 
III. FERMIONS 

A set of orthonormalized positive (negative) frequency solutions {Fa {x)} of Dirac equation in 2d MS 

I?_7/.(x)=0, V±=(i^°-^^+ij^^±7n\ , (45) 

(^0,1 _ two-dimensional Dirac matrices, 7° = /3, 7^7^ = —7^7*^ = a, a^ = /3^ = 1) is complete if 

dai^i±)(a;')Fi±)^(a;")7" = -iS^^^^x' - x"), (46) 



where the positive (+) and negative (— ) frequency Wightman functions for the fermion field are equal to 

5(±) (x) ^ 2?+ A(±' (a;), A^") (x) = A W* (a;) , (47) 

and A'^"'"^(a:) is defined in ([5]). 

Any set {Fa (x)} can be constructed of the set of matrices {S^^\x — u)} in the perfect analogy with the boson 
case. Taking into account that the Dirac operators 2?± commute with the fermion boost operator Bj}, 

[D±,BvU =0, Bv = B,c- ^iV = i [zdt + td, - ^a") , (48) 

it is easy to check that functions 

00 
i,f\x) = j dqS^^\x - u,)fM (49) 

— 00 

{uq are defined in (fTTI) ) are the eigenfunctions of this operator if 

/„(g) = e-««+^"?C, (50) 

where C is an arbitrary constant column. Substituting ()50p into Eq. (|49p we easily get with the account of definition 
(|Tfl) the following expression for the normalized boost modes 

00 

— 00 

Here the arbitrary constant column -q satisfies the relation 

?7+(l±/3)?7 = l. (52) 

It is difficult to compare representation (|51l) with the results of Refs. [14l - [l6| directly since the authors of these 
works by unknown reasons use in 2d spacetime 4x4 matrices which constitute a reducible representation of the Dirac 
matrices. The passage to their representation may be realized as follows 

a^a3=a(g)CT3, (3^(3(8)1, V ^ V^ i n) ' (^3) 



,0 
where 



-^;j). ^^(;-"0-^;^(;). (-> 



In particular, using the notations of Refs. [ij, [l| 



1 , Xa 



f = — In ■ 



2 ~x- 



V^+i-x-)^ H^UJ, $^ = ij£l^^/2(iTOu), (55) 



for the positive frequency solution in the Rindler wedge of 2d MS we obtain from Eq. ((5T|) with account of (l53l) . (|54 



i^^d{x+)d{-x^) = jyMe~ 



V 



+ e 2' 






{x+)e{-x-). 



(56) 



It is worth noting that the resuhs of Refs. PJ-UJ], see Eqs. 21.64, 21.65, 21.77 and 21.78 in [15[, differ from ^ by 
the absence of factors e^'" in front of the functions $J in ([56l) . In the framework of the method used by the authors 
of Refs. [iJ, [iSl, the latter should have appeared as a result of change of Rindler {u,v} by Cartesian coordinates 



^{t,z) 



==e^"Xu,w) 



(57) 



see, e.g., [2|. 

Besides, the modes in Refs. [M-^i3] differ from ([55]) by additional normalization factor l/-\/3. This discrepancy has 
arisen due to incorrect method of normalization used by the authors of Refs. JJrd6j]. Indeed, they have represented 
the normalization integral in 2d MS as a sum of four terms each equal to the normalization integral in one of four 
wedges of MS. However, the correct normalization procedure implies integration over some Cauchy surface in MS. 
The most convenient surface for this purpose is the surface t — 0, compare, e.g., [9|, II4I. Then only the solutions in R 
and L wedges of MS contribute to the normalization integral and the normalization factor coincides with that one in 



From now on, we will use a representation of Dirac matrices different from (j54p , 



7° = ^ = cri,7^ = icr2,7°7^ =a^ -0-3, 



(58) 

where ai, i — 1,2,3 are Pauli matrices. Use of representation ((58|) allows essentially simplify the form of solutions of 
Eq. (|^5|) . We will see below that it also is very convenient for the procedure of transition to the limit m —>■ 0. With 
this set of Dirac matrices the representation ([ST|) for the boost modes reduces to 



i^^^H^) = i 



dqe 



^im{(t^i(T) cosh q—z sinh q)—iKq ( ^ 

±6^/2 



(59) 



Using exactly the same procedure as for the scalar boost modes ((20)) we can transform Eq. ([59)) to the form 

/ f x^ ^ la ' 

V,(±)(x)^i,/^ 






I+i- 



(60) 



7 



where again w± ~ "i-\/e^™(a;+ =F ia){x- =F ia) . 

In perfect analogy with Eq. (|32) the positive frequency boost mode (jBU)) can be rewritten in the form 

i;i+\x)=0{-x+)e{-x^f^i+\x) + 0{x+)e{-x_)^^i+\x)+e{x+)e{x^fi:i+\x)+0{-x+n^ (61) 

where in P wedge 



V-Hx) = ^,d^e- 






-X- 



-x+ 



4+f K 



H 



(1) 
1/2+i 



(my^{-x+){-X-yj 



\ 



4 + 2*^ 



H 



(1) 



(62) 



in R wedge 



V^Hx) 



e 2 



V 



'X_ 



x+ 
— x_ 
x+ 



T + ryK 



l/2+« {mVi-X+)i-X-))^ 

\ 
J 



Ki/2+iK (mVx+(-x_) 



K_l/2+^K [my^X+{-X-)j 



(63) 
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in F wedge 



^^W(x) 



2V 2 



e 2 4 



a;- \ '' ^ „(2) 



^i/2+«« {m^x+x-) 



\ 



1 + o '^ 



i/ 



-l/2+^« {m^/^T^)^ 



(64) 



and finally in L wedge 



Vl+n^) = ^yfe---^- 



-a:+ 



4+^k; 



K 



\ 



\ 



l/2+iK \my/{-X^)X^ 

^-1/2+JK yTi\/{-x+)x^j 



(65) 



(66) 



A similar representation is valid for the negative frequency modes i/'k {x) as well. 
The boost modes (|59|) . (|60)) are orthonormalized by the condition 

oo 

(V'it\ Vi^)) = I dz^it)^a;)V'l±)(x) ^ <5(.. - n'), 

— OO 

constitute a complete set in the sense ([151) . and thus can serve a basis for quantization of the fcrmion field in MS, 

oo 

^{x) = j dK {b,i;i+Hx) + btyji-H^)) , (67) 

— oo 

where &« and 6^ are annihilation operators for boost fermion particles and antiparticles respectively, which obey the 
standard commutation relations. By the same reasons as in the boson case, boost quantization of the fermion field is 
unitary equivalent to quantization in the plane wave basis 



^{x) ^ f dp (ap7^(+)(x) + d+4-Hx)) , 



where 



^P "• > \^±^ep + p) 






Vp^ 



m^ 



(68) 



(69) 



It is not difficult to ascertain that the boost modes ([5^ are connected with the plane waves (|69l) through the following 
integral transformation 



00 
V;(±)(a;)= / — ^e-"?''Vi'^)(a;), Q = arcsinh^ 

J Jl-KEr, ^ m 



(70) 



The modes V'p {x) and ^pn {x) are distributions with respect to variables p and k. This means that they define 



linear functionals 



;?W(0,x)= / dp^pl^\x)g{p), 



(71) 



d^^\i),x)^ / d^i^i^\xMK), 



(72) 
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on some sets of test functions q{p), [)(«) respectively. Such functionals naturally appear when we calculate, e.g., 
matrix elements of the field operators (|67p . (155)) between the vacuum state |0m) and the states which are particle (or 
antiparticle) wave packets, 



\ls)^ J dpg{p)a+\OM), |1„) = y rf« [)W&+|Om). (73) 

— oo — oo 

Physical one-particle states must be normalized, 

CXD OO 

{1,\1,)^ I dp\g{p)\' = l, mm) = I cIkMk)]^ = 1. (74) 

— oo — oo 

Hence, g(p), ()(/«) must be square-integrable functions of p and k respectively. Then it is easy to see that matrix 
elements (|7T|) . (|72|) are square-integrable functions of x. It is clear that functionals ([71]) . ([72]) are the elements of one 
and the same space of functions. Hence, there exists one to one correspondence between the test functions g{p) and 
f)(«): 

/inq 
-j=^Q{p)dp, (75) 

— oo 

where q = arsinh(p/r7T,) is rapidity. For the sake of convenience we will confine ourselves to g(p) belonging to the 
class of continuous piecewise smooth functions descending faster than |p j"-'^ at |p | — )■ oo. This requirement guarantees 
finiteness of the mean value of the energy of the one-particle state |lg) ([73]), 



£■(1) = / dpep (Ig I a+ap I Ig) = / dp£p\Q{p)\^, £p = vV+™^, (76) 



To ascertain the properties of test functions \){k) we will use the relation ((75|) . 

First, we will illustrate by a simple example that, unlike the case of the functional (|7ip . we need know the behavior 
of test functions \){k) for the functional ([TS]) not only on the real axis but also in the complex plane. Consider 

fl(p)^ (0 ''/2u/4 ^~^^^' ">0, (77) 

(27raep^)-^'* 

which evidently satisfies the above formulated requirements to functions q{p). Using Eq. (j75p we obtain 



2a 



TT 



1/4 



(,(k) = — e-"" +'^'=. (78) 



Direct calculation of the integral ([72]) at the point x — then yields 



,...,,„,. ,/^fH£)"'-:*V,/=(Ji'Z?,*>). m 



So we see that the value of the functional ([7^ at the vertex of the light cone is determined by the values of the test 
function ()(«) at imaginary points k, = ±i/2. 

Let us now revert to Eq. (|75p . After the change of the variable of integration p — m sinh q it can be written as a 

/ TTL cosh Q 

Fourier transform of function Q{q) — \ ^{nismhq) 

V 2tt 

oo 

!,(/.)= f dqe^^'^Biq). (80) 
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i/^+^/2l^S-^/^k 



j/o + it lm/« 



:ct/2 



ReK 



b) 



-(!/2 + ic 



X-i/'2 



ReK 



-(;/2 + jf/2) -(t/2 + ie) 
FIG. 2: Contours of integration for the upper (a) and lower (b) components of 5 (f),0). 



Under our assumptions 

0((Z) = o(e-(HOkl) ^ e>0, (81) 

as |g| — ?> oo . Then, according to the Paley- Wiener theorem [25[ the test functions ()(k) (|M)) are analytic in the strip 

- 1/2 - e < Im K < 1/2 + e, e > 0. (82) 

Let us calculate the integral in Eq. (|M| by parts. Since the test functions g(g) are continuous and descend 
exponentially as \q\ — >■ oo, see Eq. (|5T|) . the first integrated term is equal to zero. Hence, the functions (}(«;) in their 
domain of analyticity ([5^ descend at least as k^^ when |k| — >■ cjo . This means that the path of integration in ([7^ 
can be shifted in the range of the strip (|82|) . 

Consider the functional (I72p at the point a; = 0. Shifting the path of integration for the upper component upward, 
and for the lower one downward, by i/2, we obtain 



;?(^Hf),o) 



27r 



dndqe 



-inq 



\){K + i/2) 
±\){K-i/2) 



: h^'i4S:-'i: 



t)(V2) 
±[)(-z/2) 



(83) 



Eq. ((83|) generalizes formula ((79)) to the class of functions [)(«;) analytic in the strip ([82)) and descending rapidly enough 
when |k| — >■ oo. 

The obtained result makes it possible to write down the fermion boost mode in the vertex of the light cone in terms 
of 5-functions of complex argument 

v-?'(o)^;f(,t-;f;^,)^ («^' 

The functionals 3^^^)(t),0) can be also written as integrals around the closed contours represented in Figs. [2^, b 

5(±)([),0)- (^ dKyji^>{OMK). (85) 



Ca,b 



/,(±) 



Indeed, the explicit form iorjp^ (0) can be easily derived from Eq. (jSO)) . Putting a;+ 
ascending series for Ki,{z) [2^, we get 

r{l/2 + iK) 

|.l/2 + iK 




^i^^(0)=lim:^| 'l'^y.r({l2-n) 



there and using the 



(86) 



yl/2-iK 



Then, taking into account that the integrand in (|55)) has simple poles at k = ±i/2 and using the Cauchy's residue 
theorem, we immediately reproduce the result (|83)) . Since the integrals taken along the vertical segments of the 
contours Ca,b are evidently equal to zero, this means, in particular, that the integrals along the upper (lower) piece 
of the contour Ca (Cf,) is also equal to zero in the limit cr — > 0. The latter statement can be easily checked by direct 
calculation. 
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^-function of complex argument was first introduced by Gel'fand and Shilov in Ref. [2Jl on the basis of Cauchy 
residue theorem, 



{Sin - ..o), fin)) = -^ / ^^^^ = fino), 
Zm J K — kq 



(87) 



where L is a contour enclosing an arbitrary complex point kq. The distribution (j87p was defined in [23'| on some 
class of entire functions /(k), the so called Z-class. Our class of test functions analytic in the strip ((82)) differs from 
Z. Therefore some properties of the Gel'fand's ^-function [23| cannot be applied to our one. However, the opposite 
statement is true: all properties of our ^-function hold valid for the Gel'fand's one. 

So, on the class of test functions analytic in the strip ([5^ . as well as for the Z-class, the following equality of 
distributions is valid 11311 



,. 1 T{l/2±iK) ^. ... 



(88) 

Let us show now that the modes ([59]) , ((60)) possess J- function singularities not only at the vertex of the light cone 
but at lines a;± = as well. Indeed, at the surface of the light cone the arguments of Macdonald functions in (|60)) are 
small, \w±\ <^ 1. Therefore using the ascending series for Macdonald functions [22] we get 



^i^H^) 



o±7rK/2^i7r/4 
2^ 



/ 



v 



/TO 



^« r(l/2 + iK) 
(a;+T»o-)i/2+i 
ril/2 -iK) 
i-x^ ± ia)^/^- 



IK.+ l 



+ i 



(f) 



ri-i/2-iK) \ 

(-a;_ ± icr)-i/2-iK 
-^K+i ri-l/2 + iK) 

(x+=F«cr)"^/^+*''/ 



(89) 



/.(±)^ 



So we see that at the surface of the light cone ipK (a;) contains (5(k — i/2) at a;+ = in the upper, and (5(k + i/2) at 
x_ = in the lower component. 

It was shown in the preceding section that in the case of a scalar field the family of boost modes does not constitute 
a complete set in MS after excluding the zero mode. This is because the scalar boost modes possess a ^-function 
singularity at the surface of the light cone, and hence the point k = gives finite contribution to physical quantities. 
Some of them, e.g., Wightman function, are determined by the spectral point n = entirely. There is a question 
whether it is possible to delete the point n = from the spectrum in the fermion case. To clarify this issue we will 
consider the integral 



■ Kq — S 



KO+S 



3^(±)([,,0) = lim / + / \dKil;i^\0)t){K)^^^'^\i),0)-\im / dKil^i^\0)t){K) . 

\ — CXD Ko+S / KO— <5 

where kq is an arbitrary real number. 

Consider the upper component of the second term in the RHS of Eq. ()90p 



(90) 



Ko+S 

^limlim f d.(^y'' '^^'ilt''\ in). 

2TTS^aa^0 J \2) crl/2+« '^ ' 

Kq-S 



(91) 



The integral in (jHT)) can be easily calculated in the limit c — > and we obtain 

«0+'5 -I/O 

(^yj r(l/2 + mo+25)t)(Ko + '^)'T 



hm y rf«[V^i^^(0)]„p()(^)=2^^2 



lini lim 



5^0 cr-i-o In a 



-tS 



(^Y r(l/2 + iKo - i5)\)ino - 5)a'' + O (-^ 



(92) 



Eq. (|92)) evidently shows that the result of calculation for S^Lp (f),0) strongly depends on the sequence of limit 
processing 5 — >■ and ct — >■ 0. If we changed the sequence of limit processing in (jOT)) . (|M)) . we would return from 
^^=^Htl.O) dSl) to S'(±)(f'>0) ([72]). Not surprisingly, Eq. ([90]) reproduces the result dH]) in this case. The accepted 
sequence which corresponds to calculation of the functional ()90p leads to a meaningless result. This means that 
pricking of any point out of the real axis k, is inadmissible in the fermion case. In other words, the set of boost modes 
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([59)1 ■ (|60)) does not constitute a complete set in MS after deleting a single arbitrary point from the spectrum. It can 
be complete only in the space which is MS without the line cone (four separate wedges of MS). In this sense our 
conclusion repeats the analogous result for boson field discussed in the preceding section. 

It is worth noting that the latter result was obtained under minimum possible restrictions as regards to the test 
functions q{p). If we restrict ourselves to g(p) belonging to the class K of smooth functions with compact support, 
as it is often done in quantum field theory, functions f)(K) would belong to the class Z [23|. In this case the proof of 
impossibility of deleting an arbitrary point kq from the spectrum is especially simple. Indeed, the Gel'fand's (5-function 
is an analytic distribution in the whole complex plain k. Thus it can be expanded in a Taylor series around any point 

K — Kq, 

-E^Ti^oT^ S'-^H^-^o), (93) 

n=0 ' ^ ^ 

and the radius of convergence of the series (|93)) is equal to infinity [23]. Here (5'"^ (k — kq) denotes the n-th derivative 
of the (5- function at the point k — kq. As a result, 



dK 5 ("k T ^) ()(«) = ^ ^ Uo T ^) Irni J + y dnS'^^^ {k - Ko)f)(«) = , (94) 




K-O + i ' 



and hence, on the class Z of test functions, deleting of an arbitrary spectral point kq from the spectrum leads to 
vanishing of functionals 5^(*)((), 0) ([85]) at the vertex of the light cone, and its finite variation at other points of the 
cone surface. However, since the Gel'fand's (5-function is analytic on the whole plane of complex k, it can be expanded 
in a Taylor series around another real point k — ki, ki ^ kq- If we calculate integral (j94p using such representation 
of the (5-function, we will see that deleting of the spectral point kq from the spectrum does not influence the result 
of integration. Thus, since the value of the matrix element ()85p depends on the way of calculation, we conclude that 
the operation of exclusion of an arbitrary point from the spectrum is meaningless in agreement with the previous 
consideration, see discussion of Eq. (|90p in the preceding paragraph. 

As it was shown in the preceding section, Wightman functions A^*' [x) of a massive scalar field are determined 

by the zero boost modes ^q (x), see Eq. p5|) . This is a consequence of translational invariance of A^^^(a;) and the 
presence of Dirac J-function singularity of ^L (x) at the light cone. A similar result is valid also for the Wightman 
functions of a massive fermion field. Indeed, due to the property of translational invariance we can write, compare 
dUl), 

GO 

5(±)(x',x") = i / dKVi^)(a;)V'i^)'(0)7°, x = x' - x" . (95) 



Using now Eq. (|84p we get the following result for matrix elements •S'^a (x) of Wightman function (IMj) . 

S'^^\x)^±^^i^fl{x), S^M\^)^^\I^^%S^), " = 1,2, (96) 

where V'k;,q(x) are the a-components of the boost modes ([59|) . ((60|). Thereby, we see that the matrix elements of 
S'^^^{x) are determined by only two "spectral points" k = ±i/2. Exactly as in the scalar case, it is easy to ascertain 
that this result holds valid for the smeared Wightman functions as well. 

Using Eq. (|60l) one can express the Wightman functions (l96l) in the form, compare (j44|), 



^(±)(x)=^f 

ZTT 



^ ±Ko(-±) ('-±^''' K.i^n^^ 

\x^ =p l(J J 



(97) 



where w± is the same as in Eq. ([m]). 

Consider now the case of a massless fermion field first studied in Ref. [l^ . The Wightman function of the massless 
field Sq (x) can be obtained easily by passage to the limit m ^ in Eq. (1^7]) 

/ ~ 1 



1 ^ 

S^,^\x)^ lim^(±)(x)=±- 



X- =F la 



^+ ^ *'^ I , a^O. (98) 
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Just as in the cases considered earlier, one can use 5'g (a; — Uq) to obtain any complete set of orthonormalized 
solutions of Dirac equation. In particular, for the boost modes we have, 

oo 

xtHx)= I dqS\^\x-aq)e-^-'i-i^^C, (99) 



where C = { j is an arbitrary two-component column, compare (|49|. (l50l) . We see that at to = the Wightman 

function (|98p is antidiagonal. Owing to this property the upper and the lower components of the boost mode ((99|) 
become independent: 

Xl±^ {x) = £ixit' (2^)e+ + 52x1^^ (2^)e- , (100) 

where e± are the eigenvectors of the Pauli matrix (T3, a^e± = ±e+ . This is because there appears a new conservation 
law for the massless case, conservation of chirality, see, e.g., [2^|. Thereby solutions of the Dirac equation (|45l) are 
labeled by the additional quantum number r = ± , chirality. The coefficients ci and £2 in (llOOp are independent and 
can be determined by the normalization condition for solutions x];.^ (x)e+ and xj,_ (a;)e_ separately. Note that in the 
massive case, since chirality is not conserved, the components of the column C (|50p cannot be determined separately. 



In that case their combination forms a monomial factor of the solution V'k {x) (|51|) . 
For the normalized positive frequency functions Xkt we have 

where the phase factors were chosen for the sake of convenience. It is worth noting that the negative frequency 
functions Xkt come out of p01[) by complex conjugation and the change k — > —k. 
Plane waves are also labeled by the quantum number t in the massless case, 

^(t'(•^)-4=eT'^^-^ 0<p<^. (102) 

They are linked to the boost modes Xkt (a;)eT through the Mellin transform 

CXj 

xW(x.)=|-^/-V^^Hx.). (103) 



This means that, if we assume that distributions (ppr {xr) are defined on the same class of test functions q(p) as in 
the massive case, the distributions Xht {xt) will be defined on the class of test functions \){k) analytic in the strip 



and descending at |k| — > 00 in this strip. 
Then taking into account Eq. (IM|) we conclude that boost modes of a massless fermion field at the surface of the 
light cone are (5- functions of a complex argument, compare |13l |. 



x'■^,Kxr = Q)^5\^^^-T-y (104) 

Now we will show how the modes (jlOip can be obtained from Eq. (j60p by limit processing to — J- 0. Using the 



ascending series [22| for functions Ki/2±iK{w±) we get 

^^^Hx) = (^)-^\it^e,±(!^)^M-Je_, (105) 

i.e., representation (|100p with coefficients ci and £2 containing singular at m ^- phase factors. However these factors 
do not influence the normalization constants, have no impact on any physical quantities and hence can be omitted. 
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IV. CONCLUDING REMARKS 

We have shown that Wightman function of a free quantum field generates any complete set of solutions of relativistic 
wave equations. Using this approach we have constructed the complete sets of solutions to KFG and Dirac equations 
consisting of eigenfunctions of the generator of Lorentz rotations (boost operator). 

Boost modes are used as a basis for field quantization very rarely. Till now they were exploited only for analysis of 
the so-called "Unruh effect" [9, 12 1 and at attempts to quantize a charged massive scalar field in the presence of an 
external constant electric field 



u 



28| . However, there are many problems, especially in the quantum field theory in a 
curved space, where the boost symmetry may appear to be the only symmetry for the quantum field, and thus using 
it for separation of variables in a classical field equation is the only instrument to find solutions for such equations. 

The specific feature of the boost modes is that, taken at the surface of the light cone, they as functions of the boost 
quantum number k possess strong singularities. Certainly, the physical reason for these singularities is the singularity 
of Lorentz transformations at v = c. 

For the case of a scalar field this is a Dirac (5-function singularity ^^(2^)1 _r, ~ S{k) 12]. This leads to a special role 
of the zero boost mode: its exclusion from the set of the boost modes makes the latter incomplete. The exceptional 
role of the zero boost mode becomes quite clear if we recall that the zero value of k means that ^0(2;) is a Lorentz 
invariant positive frequency solution of KFG equation, i.e. the Wightman function for the quantum field coinciding 
with the positive frequency part of the commutator of two scalar field operators. Hence the exclusion of the zero 
boost mode results in a "quantum" theory with commuting field operators. 

The singularities of the fermion boost modes are even stronger. It is shown in the present paper that at the surface 
of the light cone they possess 5-function of a complex argument ^KQ (2:) _(, ~ (5(K=pi/2). The ^-function of a complex 
argument was first introduced by Gel'fand and Shilov in Ref. [23|, and was defined on the class Z of entire functions. 
Our (5-function is defined on the class of test functions analytic in the strip ((82)) . Actually, the width of the strip 
is determined by physical requirements. Our choice provides square integrability of one-particle wave packets and 
finiteness of their energy. If we require finiteness of the squared energy we should narrow the class of test functions 
and extend the width of the strip to 

-3/2 - e < ImK < 3/2 -I- e , e>0. 

Further toughening of requirements to physical states will lead to subsequent extension of the width of the strip 
of analyticity of the test functions. Therefore it is reasonable to generalize the concept of (5-function of a complex 
argument and introduce the distribution (5s(k) 

iSs{n-Ko)J{K))=fiKo), (106) 

defined on the class Zs of test functions /(k) analytic in the strip 

- s - e < ImK < s + e, s > 0, e > , (107) 

and kq belong to this strip. In such notation the (5-function introduced in Sec. HI will look as Si/2{k), the Gel'fand 
(5-function as Sooin) and the standard Dirac 5-function as 6o{k). 

It is worth noting that for (5-function defined on the class Zs the following integral representation is valid 

00 

(5^(zy) = -L f dqe-""^, i^ = KTis, s > 0, (108) 

2tt J 

—00 

which at s = is a straightforward generalization of the standard representation for the Dirac (5-function. Besides, 
the following relations are valid, 

'5.(^) = ^4!^(0)--^i/L''(0) = -if,.(0), ^^KTis, s>0, (109) 

compare Ref. ^2^] for the case s = 0. 

The presence of a 5-function of a complex argument in a boost mode at the surface of the light cone does not allow 
to exclude any point from the boost spectrum. We have shown this for the case s = 1/2, see Eqs. (P(I1) - ([M1) . and this 
statement proving does not change for an arbitrary value of s. Another way to prove this statement was used for the 
distribution 5oo(k) based on its analyticity in the whole complex plane. The analogous proof could be given for the 
5s-function as well. In this case, due to the finite value of radius of convergence R = s + e, the expansion (|93l) can be 
applied only for real kq = bi , |&i| < e. As a next step the Dirac 5-function 5o(k — 61) and every its derivative can be 
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represented as a Taylor series of the type (l93l) centered at the point 62 on the real axis, I62 — &i| < s + e. So, Ss{k — is) 
will be expanded in a 2-multiple series of the Dirac 6o{k — 62)-function and its derivatives. After a finite number of 
steps N we can reach an arbitrary spectral point kq and hence obtain a representation oi 5s{k — is) in the form of a 
A'^-multiple series of the Dirac ^o('*~ Ko)-function and its derivatives. It is clear that the point kq cannot be excluded 
from the spectrum then. 

The latter reasoning clearly explains why the scalar case, when we have the only one distinguished point k = 
which cannot be excluded from the spectrum, differs drastically from the fermion case. Indeed, the strip of analyticity 
(I107P for Ss{k) degenerates into the real axis of the k complex plane at s = 0, so that the radius of convergence for 
the corresponding Taylor series becomes equal to zero, or in other words the Dirac (5o('^)-function is not an analytic 
distribution. Thus the procedure discussed above cannot be realized. 

To conclude, it is worth emphasizing that we have shown explicitly that smearing of boost modes, or Wightman 
functions does not change our results, see also Ref. |3(j |. 

V. ACKNOWLEDGEMENTS 

This work was supported by Russian Fund for Basic Research, the Ministry of Science and Education of the Russian 
Federation and the Russian Federal Program " Scientific and scientific-pedagogical personnel of innovative Russia" . 



[1] 
[2] 
[3] 
[4] 



[5] 
[6] 
[7 

[s; 
[lo; 

[11 

[12 

[13 
[14 
[15 
[16 
[17 
[18 
[19 

[20 

[21 
[22' 
[23^ 
[24 
[25| 

[26 

[27 
[28' 
[29' 
[30 



N.N. Bogoliubov, D.V. Shirkov, Introduction to the theory of quantized fields, NY.: John Wiley (1980). 

A.I. Akhiezer, V.B. BeTestetskii, Quantum electrodynamics, NY: Interscience Publishers (1965). 

A. I. Nikishov, Zh. Eksp. Teor. Fiz. 57, 1210 (1969) [Sov. Phys. JETP, 30, 660 (1969)]. 

N. B. Narozhny, A. I. Nikishov, Theor. Math. Phys. 26, 9 (1976); Issues in Intense-Field Quantum Electrodynamics: 

Proceedings of the Lebedev Physics Institute of the Academy of Sciences of the USSR: v. 168, ed. by V.L. Ginzburg (Nova 

Science Publishers, New York, 1989) p. 226. 

D.M. Volkov, Z. Phys., 94, 250 (1935). 

P. J. Redmond, J. Math. Phys., 6, 1163 (1965). 

D. Grumiller, W. Kummer, D. V. Vassilevich, Phys.Rept. 369, 327 (2002). 

S.V. Christensen, S.A. Fulling, Phys. Rev. D 15, 2088 (1977). 

W.G. Unruh, Phys. Rev. D 14, 870 (1976). 

A.I. Nikishov, V.I. Ritus, Zh. Eksp. Teor. Fiz. 94, 31 (1988) [Sov. Phys. JETP, 67, 1313 (1988)]. 

U. Gerlach, Phys. Rev. D 38, 514 (1988). 

N.B. Narozhny, A.M. Fedotov, B.M. Karnakov, V.D. Mur, and V.A. Belinskii, Phys. Rev. D 65, 025004 (2002). 

A. M. Fedotov, N. B. Narozhny, V. D. Mur and E. G. Gelfer, JETP Letters 89, 385 (2009). 

M. Soffel, B. Miiller, W. Greiner, Phys. Rev. D 22, 1935 (1980). 

W. Greiner, B. Miiller, J. Rafelski, Quantum Electrodynamics of Strong Fields, (Springer- Verlag, 1985) 

D. McMahon, P. M. Alsing, and P. Embid, arXiv:gr-qc/0601010 

R. Jost, The General Theory of Quantized Fields, (AMS, Providense, 1965). 

R.F. Streater, A.S. Wightman, PCT, Spin and Statistics and All That, (Princeton University Pres, Princeton, 2000). 

N.N. Bogoliubov, A. A. Logunov, I.T. Todorov, Introduction to axiomatic quantum fields theory, (W.A. Benjamin, Inc., 

New York, 1975). 

A.S. WightmsLTi, Introduction to Some Aspects of the Relativistic Dynamics of Quantized Fields, (Princeton University Press, 

Princeton, 1964). 

L.C.B. Crispino, A. Higuchi, G.E.A. Matsas, Rev. of Mod. Phys. 80, 788 (2008). 

I. S. Gradstein and I. M. Ryzhik, Table of Integrals, Series, and Products, (Academic Press, New York, 1965). 

I.M. Gel'fand, G.E. Shilov, Generalized functions, Volume 1, (Academic Press, New York, 1964). 

V.A. Fock, Zeit. Phys. 57, 261 (1929). 

R.E.A.C. Paley and N. Wiener, Fourier Transforms in the Gomplex Domain, (American Mathematical Society, New York, 

1934). 

N. V. Krasnikov , V. A. Matveev , V. A. Rubakov, A. N. Tavkhelidze, V. F. Tokarev, Theor. Math. Phys. 45, 1048 (1980) 

C. Gabriel and P. Spindel, Annals Phys. 284, 263 (2000). 

N.B. Narozhny, V.D. Mur and A.M. Fedotov, Phys. Lett. A 315, 169 (2003). 

D.V. LVov, A.S. Shelepin, L.A. Shelepin, Yad. Fiz. 57, 1147 (1994). 

N.B. Narozhny, A.M. Fedotov, B.M. Karnakov, V.D. Mur, and V.A. Belinskii, Phys. Rev. D 70, 048702 (2004). 



